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^SJ , Using the covariant formalism, we derive the equations of motion for adiabatic and entropy 

Oh. perturbations at third order in perturbation theory for cosmofogical models involving two 

• scalar fields. We use these equations to calculate the trispectrum of ekpyrotic and cyclic 

models in which the density perturbations are generated via the entropic mechanism. In these 
models, the conversion of entropy into curvature perturbations occurs just before the big 
bang, either during the ekpyrotic phase or during the subsequent kinetic energy dominated 
^ • phase. In both cases, we find that the non-linearity parameters Jnl and gj^^ combine to 

pH ' leave a very distinct observational imprint. 
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I. 



INTRODUCTION AND SUMMARY OF RESULTS 



The cosmic microwave background radiation (CMB) provides us with a staggering amount of 
information about the early universe. In fact, its exploration has only recently begun in detail 
with the WMAP satellite, but with the successful launch of the Planck satellite, we can hope to 
get yet more precise maps of the microwave sky in the near future. At the same time, surveys 
of the large-scale structure of the universe are becoming increasingly detailed, to the point where 
they can rival the precision of the CMB maps, albeit with more complicated, and less understood, 
physics involved. The analysis of the statistics of density fluctuations is thus becoming a strong 
discriminator of competing models of the early universe. 

An increasingly important aspect is non-gaussianity [1], i.e. the deviation from perfect gaussian 
statistics. Indeed, there are many models of the early universe that, at the linear level, give 
degenerate predictions {e.g. the large class of inflationary models [2] and ekpyrosis [3]) - hence 
the importance of studying higher correlation functions. As we go to higher n-point functions, the 
degeneracy typically tends to be broken, as has been demonstrated by the study of the bispectrum 
in recent years. Clearly, the larger the number of distinct observations that a cosmological model 
is in agreement with, the more plausible it becomes, and this constitutes the primary reason for 
extending the analysis to the trispectrum, i.e. the Fourier transform of the 4-point function. 

Even though we are looking at what appear to be increasingly small details of the CMB maps, 
the implications can be far-reaching, since difl^erent models can give rise to vastly different pictures 
of the universe in terms of its long-time evolution and its appearance on scales larger than the 
current horizon - in this vein one should for example compare eternal inflation [4] and the phoenix 
universe [5]. Hence it is of crucial importance to understand the predictions of the various models, 
ideally before the measurements are made, so that these measurements can be interpreted within 
a suitably developed theoretical framework. 

In this paper, we extend to third order in perturbation theory the derivation of the equations of 
motion needed to calculate the non-gaussianity of two-scalar-field cosmological models. For this, 
we use the covariant formalism, as it provides a conceptually unambiguous path to the definition 
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of higher-order perturbations and the derivation of their equations of motion (for some of our 
examples, the background evolution is complicated, so that we cannot easily use the SN formalism). 
The resulting equations are generally applicable, and we have used them here to calculate the 
trispectrum of perturbations generated by ckpyrotic and cyclic cosmological models. For certain 
inflationary models, the trispectrum has already been calculated using the 6N formalism [6-8], 
and, in future work, it would certainly be of interest to use the formalism developed here in order 
to check and extend these results to cases where the dN formalism is less readily applicable. 

Ekpyrotic models constitute a case of particular interest regarding non-gaussianities, since the 
cosmological perturbations are generated in the presence of a steep potential. Hence, we are 
automatically in a regime where we expect significant self-interactions of the scalar fields, and 
hence significantly non-gaussian distributions of temperature fluctuations. This is in stark contrast 
to simple inflationary models [9, 10], where the flatness of the potential implies that the inflaton 
is practically massless, with negligible sclf-intcractions. In more complicated inflationary models, 
such as the curvaton model [11, 12], significant non-gaussianities can be developed due to the 
interactions of various scalar fields, but with the drawback that the range of possible values for the 
non-linearity parameters is difficult to pin down [65]. The beauty of the ekpyrotic models is that, 
despite the fact that the potential has not been derived microphysically, and that we therefore must 
use a parameterized form for the potential, the results are nevertheless surprisingly predictive. 

In this paper, we restrict our attention to those models in which the curvature perturbations 
are generated via the cntropic mechanism [13], as this is the most robust mechanism known to 
date. In this mechanism, nearly scale-invariant entropy perturbations arc generated first, and are 
subsequently converted into curvature perturbations via a bending of the trajectory in scalar field 
space just before the big bang. Such a bending of the trajectory is generically expected in various 
contexts [14-16]. Broadly speaking, there are two limiting cases that are of special interest: the 
first is where the bending occurs after the ekpyrotic phase, during the approach to the big crunch. 
In this case, which is the natural case in cyclic models, the kinetic energy of the scalar fields is the 
dominant contribution to the total energy density while the conversion takes place, and hence we 
call this case "kinetic conversion". The other limiting case is where the conversion occurs during 
the ekpyrotic phase, and we call this "ekpyrotic conversion" . 

Postponing all details to sections IV and V, we simply state our main results here. We define 
the local non-linearity parameters /^vl and g^L via an expansion of the curvature perturbation ( 
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in terms of its linear, gaussian part (l, 

C = Cl + ^/jvlC! + ^9nlCI (1) 

During the ekpyrotic phase, we adopt the following parametrization of the potential: 

Ve, = -Vbe^-[1 + es' + ^e'/'s' + '^eh' + •••], (2) 

where we expect K3, K4 ~ C(l) and where e ~ O(IO^) is related to the ekpyrotic equation of state 
via e = 3(1 + Wg,k)/2. a denotes the adiabatic direction, i.e. the direction tangent to the scalar 
field space trajectory, while s denotes the "entropy" direction, i.e. the direction perpendicular 
to the background trajectory. During the ekpyrotic phase, we are able to solve for the entropy 

perturbation, with the result that 

with the linear, gaussian part Ssl being inversely proportional to time t. 

Then, for kinetic conversions lasting on the order of one e-fold of contraction of the scale factor, 
we find the following approximate fitting formula for the third order non-linearity parameter g^L 
(and we include /atl for completeness [17, 18]) : 

3 

Inl ~ 2'^3A/e + 5 (4) 
100,(^ + 1-?). (5) 

These simple formulae accurately reflect the scaling with e, K3, K4 of the non-linearity parameters. 
But crucially, note that when K3 and K4 are small, Qnl always negative. Hence, even though /jv£, 
is small in that case, qnl is negative and typically of order a few thousand, so that any accidental 
degeneracy at the level of /jvl between simple inflationary and cyclic models is very likely to be 
broken at the level of the trispectrum. More generally, unless I/atlI turns out to be quite large, one 
would typically expect gNL to be negative, as obtaining a positive gNL in that case would require 
unnaturally large values of K4. 

For ekpyrotic conversions, it is easier to perform the analysis in terms of the orthogonal fields 
01 and (p2, with (j)2 being the direction of the trajectory after the bending (0i and (p2 are related 
to a and s by a rotation in field space). Then, if the potential is expressed as 

Kfc = -V^ie"'^*^^ - ^2e-^^'^^ (6) 
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we find the following approximate formulae for the non-linearity parameters (where Jnl was first 
calculated in [19]): 

5 

Jnl (7) 

9NL ~ ^cf. (8) 

Here /jvl is always negative while ^atl is always positive, and both scale rather fast with increasing 
ci- Since we must have ci > 8 in order not to conflict with spectral index measurements, we can 
expect future observations to be unambiguous in their support or refutation of these models, and 
current limits on /jvl are already putting some strain on this particular mode of conversion. 

The plan of this paper is as follows: we start with a short summary of the covariant formalism, 
and of how it has been used to derive the equations of cosmological perturbations up to second 
order in perturbation theory. Then we extend these results to third order, and derive the equations 
of motion relevant for calculating the trispectrum of cosmological perturbations. In section IV, 
we define the non-linearity parameters that are typically used to quote observational results, and 
that we wish to calculate. Then, in section V wc apply our equations to ekpyrotic models, and 
numerically solve them for the cases of kinetic and ekpyrotic conversion. We end with a discussion 
of our results. The appendix contains a collection of useful formulae. 



II. COVARIANT FORMALISM AND PERTURBATION THEORY UP TO SECOND 

ORDER 

Motivated by the calculation of the primordial trispectrum in ekpyrotic models, we study the 
cosmological fluctuations of a system of two scalar fields at third order in perturbation theory. We 
adopt the 1-1-3 covariant formalism developed by Langlois and Vernizzi [20-23], along the lines of 
earlier works by Ellis and Bruni [24] and Hawking [25] (see also [26] ) . Applications of this formalism 
to a wide variety of cosmological contexts can be found in [27]. It is particularly well suited to 
analyze fluctuations of the system we consider, as it provides a geometric and, most importantly, a 
fully non-perturbative description of spacetimes that are distorted away from the standard FLRW 
form. 

The idea is to define perturbations via spatial projections of various scalar quantities, such as 
the energy density or the integrated expansion along a time-like curve. These spatial projections 
(covectors) vanish for an exact homogeneous and isotropic FLRW spacetime, and hence they nat- 
urally play the role of 'perturbations', even though they can be defined non-perturbatively. The 
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formalism allows one to derive simple, all-orders, equations of motion for these covectors, which, 
after choosing a convenient coordinate basis, can be expanded up to the desired order in pertur- 
bation theory. Also, the formalism leads to a clear and natural decomposition of fluctuations into 
adiabatic and entropy components, again in a fully non-perturbative way. 

In this section, we give an overview of the covariant formalism and its application to the study 
of a system of two scalar fields, up to second order in perturbation theory, as done in [23] (see also 
[28] for an extension to two-field models with non-standard kinetic terms of the form studied in 
[29, 30]). This provides the necessary groundwork for our investigation of third-order perturbation 
theory in the following section. We will also focus on the large scale limit, i.ewe only consider the 
leading-order terms in an expansion in spatial gradients. In this respect, our formalism is similar 
to the one developed in [31]. 

Our analysis is quite general: we do not specify an underlying potential and do not use approx- 
imations other than the large scale limit. In particular, although we will apply it to a particular 
scenario of a contracting universe, it is equally well suited to the study of cosmological perturbations 
generated during a period of inflation. 

A. Covariant formalism 

A key role in the covariant splitting of space and time in the 1-1-3 covariant formalism is played 
by the a priori arbitrary unit timelike vector u"' = dx"'/dT {uau"' = — 1), defining a congruence of 
cosmological observers. The spatial projection tensor orthogonal to the four- velocity is provided 

by 

hab = gab + UaUb, {h\hl'^ = h\, hSb = ^)- (9) 

where gab is the cosmological metric. To describe the time evolution, the covariant definition of 
the time derivative is then the Lie derivative with respect to u", defined for a generic covector Ya 
by (see e.g. [32]) 

Ya = CuYa = u^VbYa + y^,V„u^ (10) 
and denoted by a dot. For scalar quantities, one simply has 



/ = u^'Vbf. 



(11) 
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To describe perturbations in the covariant approach, we consider the projected covariant derivative 
orthogonal to the four- velocity u", denoted by Da- For a generic tensor, the definition is 

DaT,^:- ^ ...hf... vdTjr. (12) 

In particular, when focussing on a scalar quantity /, this reduces to 

Daf = h^y^f = daf + uJ. (13) 

Geometrical properties of the congruence are defined via the decomposition 

VfeMa = (Jab + UJab + ^^Kb - daUb, (14) 

with the (symmetric) shear tensor and the (antisymmetric) vorticity tensor coab', the volume 
expansion, 6, is defined by 

e = VaU", (15) 

while 

o" = uVbu" (16) 

is the acceleration vector. A useful related quantity is the integrated volume expansion along u"-, 
defined as 



\jdTQ (e = 3a). (17) 



a = 

3 

Since 9/3 corresponds to the local Hubble parameter, one sees that the quantity a can be in- 
terpreted as the number of e-folds measured along the world-line of a cosmological observer with 
four- velocity n". Finally, it is always possible to decompose the total energy momentum tensor of 
a system as 

Tab = {P + P)UaUb + QaUb + UaQb + GabP + T^ab, (18) 

where p, P, Qa and -Kab are respectively the energy density, pressure, momentum and anisotropic 
stress tensor measured in the frame defined by tx". 

The above quantities enter into the definition of the generalized curvature perturbation on 
uniform energy density hypersurfaces (see [20, 21]) 

di 

Ca = Daa - -DaP. (19) 
P 
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is the key quantity we want to track the evolution of. In the large scale limit (for which we are 
using the notation ~), it satisfies the simple first-order evolution equation 

Ca ~ 3^ra, (20) 

where 

r„ = DaP - jDaP (21) 

is the nonlinear nonadiabatic pressure perturbation. For a barotropic fluid, Fa = and (a is 
conserved on large scales. The relation (20) can then be seen as a generalization of the familiar 
conservation law for the linear curvature perturbation on uniform energy hypersurfaces [33, 34]. 
Finally, note a useful property of (a- in its definition Eq. (19) one can replace the spatial gradients 
Da by partial or covariant derivatives, 

Ca = Vaa-jVaP. (22) 

The same property of course applies to similar combinations of spatial gradients, such as F^. 

B. Two scalar fields 

In this subsection, we specify the above formalism to the situation where two scalar fields 
minimally coupled to gravity fill the universe. The corresponding Lagrangian density is 

jC = -\dact>id''<t>' - V{ct>K), (23) 

where V is the potential. Wc assume, for simplicity, canonical kinetic terms. Here and in the 
following, summation over the field indices (/, J = 1,2) will be implicit. The energy-momentum 
tensor derived from this Lagrangian reads 

Tab = da4>idb4>' - \gah (^c^/^^^' + 2V) . (24) 
Starting from the energy-momentum tensor (24) one finds 

p = TabU^'u^ = ^ (j>i4>' + Da^'D''<t>i) + y, (25) 
P = ^h-'^Tabh', = ^ - ^DacPiD-cp'^ - V, (26) 

Qa = -U^chl = -4>lDa(t>', (27) 

7r„6 = KTcdhi - Phab = DahDb^' - habD,ct>iD''ct>' . (28) 
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In the two-field case, it is convenient to introduce a particular basis in field space in which 
various field dependent quantities are decomposed into adiabatic and entropy components. In the 
linear theory, this decomposition was first introduced in [35] for two fields. For the multi-field case, 
it is discussed in [29, 36] in the linear theory and in [31] in the nonlinear context. In our case, 
the corresponding basis consists of a unit vector defined in the direction of the velocity of the 
two fields, and thus tangent to the trajectory in field space, and of a unit vector el defined along 
the direction orthogonal to it (with respect to the field space metric Gjj, taken here to be trivial, 
Gij = Sij, with which the field indices are raised and lowered), namely 

ei^ ^ ei= ^ (-02, (29) 

To make the notation shorter it is convenient to introduce the angle 9 defined by 

a • a ^"^ 

COS 9 = — -j^=^= , sm ^ = — 1^=^= , (30) 

V^!+^ ^l4>l + 4>l 

so that 

= (cos^,sin^) , = (— sin6',cos^) . (31) 

This angle, in contrast with the linear theory case where it is a background quantity that depends 
only on time, is here an inhomogeneous quantity which depends on time and space. By taking the 
time derivative of the basis vectors and e^, we get 

4 = 9el ei = -9ei. (32) 

It is also convenient to introduce the formal notation 

a= VG7j<^V = -v/</'? + 0i, (33) 

such that, for instance, 

ei = t (34) 

Notice that generically, the quantity a is not the derivative along of a scalar field a ; we just 
use it as short-hand notation. 

Making use of the basis (29), one can then introduce two linear combinations of the scalar field 
gradients and thus define two covectors, respectively denoted by aa and Sa, as 

aa = eaiVact)^ (35) 
Sa = esiVa^^- (36) 
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These two covectors are called the adiabatic and entropy covectors, respectively, by analogy with 
the similar definitions in the linear context [29]. 

Although one can derive a second-order (in time) equation for the adiabatic covector (see [23] ) , 
on large scales it is simpler to resort to the evolution equation for (^a- In that purpose, it is useful 
to introduce the covector 

^a = Dap- ^a^, (37) 
a 

which can be interpreted as a covariant generalization of the comoving energy density perturbation. 
Indeed, it has been shown in [23] that if the shear can be neglected on large scales, Ca is also 
negligible, in particular in Eqs. (38) and (39) below. On large scales, the shear rapidly decreases 
in an expanding universe while it gets far less blueshifted than the dominant energy contribution 
in ekpyrotic models; therefore, one can safely neglect it in both cases. However, to accurately 
describe the evolution of perturbations in a local way on large scales, we will need to keep track of 
it, as discussed below. 

Then, in our two-field system, Eq. (20) reduces to [23] 

2 26 

Ca^^V^ea + ^VsSa. (38) 

while the entropic equation can be expressed as 

Sa + QSa + [Vss + 3^^J S„ « -'^^^a- (39) 

In the above equations, = ^i-Vj, Vg = ^l^j, Vgs = ^l^iViJ, etc. Recall that we are using 
the notation that denotes an equality on large scales. Eqs. (38) and (39) are our two master 
equations. They are exact on large scales and despite their simple appearance, they contain all the 
nonlinearities of the Einstein-scalar field equations on large scales. In the following subsection, we 
will review how the introduction of a coordinate system enables one to straightforwardly derive 
from them the linearized and second-order perturbation equations, while we derive new results 
at third order in the next section. Although the route to pursue is straightforward, intermediate 
calculations can become long. Therefore, for convenience, we have collected in the appendix various 
background as well as first- and second-order expressions that will be used in the rest of this paper. 

C. Perturbation theory up to second order 

To relate the covariant approach with the more familiar coordinate based formalism (for a recent 
review of cosmological perturbations in this context, see [37]), one introduces generic coordinates 
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x'^ = {t, x^} to describe an almost-FLRW spacetime. A prime will then denote a partial derivative 
with respect to the cosmic time t, i.e. ' = djdt, since the dot has been reserved to denote the Lie 
derivative with respect to n". Also, we expand fields as follows: 

X{t, x') =X{t)+ (5X(i) (i, x') + {t, x') + (5X(3) (i, x'), (40) 

i.e. without factorial factors, where a quantity with an over-bar is evaluated on the background. 
Here, 5X^^^ refers to a quantity that solves the linearized equations of motion, (JX(^) the quadratic 
equations, and so on. Occasionally, and only when we feel that the meaning is unambiguous, we 
drop the superscript for perturbations at linear order. Also, in the following, u"- is chosen such 
that Ui = 0. The remaining component uq is then determined solely in terms of metric quantities. 

By expanding Eqs. (35) and (36) up to second order, one finds, for cTj and si respectively, 

5(Ti = diSa, Sa = Efji^cj)^ (41) 
5si = dids, 6s = CsiScf)^ , (42) 

at linear order and 

,5(7f ) = diSa^^^ + ^SadiSs - ^Vi, (43) 
a a 

Ssf = dM^^ + ^-^d^5s', (44) 
at second order. We have used the definitions 

5(7^2) = e^i5^^^'^'^ + -^,5s5s', (45) 

2(7 

and introduced the spatial vector 

Vi = ^{5sdi5s' -Ss'diSs), (47) 

which we will discuss shortly. 

The definitions of the entropy perturbations at first and second order are chosen such that these 
quantities are gauge-invariant. Given a vector the gauge transformation it generates is defined 
by the transformation law of tensors [37] (whose coordinate functions are a particular case) 

f^e^«T. (48) 
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With the perturbative expansion ^ = 5^„^C(n)) the first and second-order perturbations of a 
tensor T are then found to transform as [38] 

STi^) ^ <5T(i) + C^^^^ T W , (5T(2) ^ (5T(2) + C^^^^ t(o) + + C^^^^ <5T(i) . (49) 

Using these relations, it is easy to verify tliat Ss^^^'^"^^ are indeed gauge-invariant on large scales. 

The adiabatic perturbations are gauge-variant quantities, and hence, their definition can be 
chosen for convenience. In fact, it turns out that an inspection of the momentum density 

Qi = -&ai (50) 

typically provides a clue as to what a convenient definition for Sa ought to be. A straightforward 
expansion yields 

Sqi = -di {a'6a) (51) 



= -di ( a'Sa^^^ + l^Sa^ + e'SaSs] - ^SeOiSa + Vi, 
\ 2 a' J 



(52) 



and hence we can see that the momentum density perturbations almost automatically vanish when 
the adiabatic perturbations do. In other words, setting the adiabatic perturbations to zero (as 
a gauge choice) almost automatically corresponds to adopting comoving gauge. The fly in the 
ointment is the spatial vector Vi. However, as can be seen from its definition above, Vi vanishes 

when 6s' and 5s have the same spatial dependence, i.e. when Ss' = f{t)5s. This condition is 
fulfilled for super-Hubble modes in inflationary and also in ekpyrotic models, as such relative 
spatial gradients are heavily suppressed on large scales in both cases. Thus we can define an 
approximate comoving gauge at second order, which coincides with setting 6a^^^ = Sa^"^^ = 0. 

Expanding now Eq. (37), one finds 

Sei = di6e, 6e = dp — zr^da, (53) 

a 

where Se is the first-order comoving energy density perturbation, while at second order, 

<5ef = 5,<5e(2) + ^-^diSe^'^' -3HVi, (54) 

with (5e(^) defined by 

Sei^) ^ 5p('^ - ^Sa^'^ - ^ L^^)' + ^(l\6a + h'5s 
a' a' \_ 2 \ cr' / a' 



(55) 



representing the energy density in the approximate comoving gauge on large scales. 
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As for the curvature perturbation on uniform energy density hyper surf aces, it reads 

SCi = diC, C = Sa- j5p. (56) 



= ^.C^^) + |5,C(^)' (57) 



;(2) ^ f.A'^) ^ 
with 

n ... An r 1 / ff \ ' 1 

(58) 



Using the transformations (49), one can readily verify that ^e^^-*'*-^^ and ^(^)'(2) are gauge-invariant 
quantities on large scales. 

After these identifications of the desired gauge-invariant quantities, the expansion of the master 
equations Eqs. (38) and (39) on large scales gives (neglecting the first- and second-order comoving 
energy density 5e and (5e(^^ at the end of the calculation) 

2 1-f 

C'^—^0'5s, (59) 
6s" + 3H6s' + {Vss + S9'^)Ss « (60) 



at first order, and 



,(61) 



^ ' ' a' a' \ a' 2 J 



1 - 6' - 6'^ \ 

^Vsss - ^-Vss - 9— Ss"" (62) 

at second order. It is worth commenting that if one neglects the covector €„ in Eqs. (38) and (39) 
from the beginning, one can write the second order equations of motion derived from them as a 
total gradient, and hence in a local form as in Eqs. (61) and (62), only if one assumes throughout 
that Vi = 0. If ta is kept however, the equations of motion become local, and the two limits of the 
equations agree if the scalar quantities (5e and 8e^^ are neglected, as is appropriate. 

These equations form a closed system: on large scales, the entropy perturbation evolves in- 
dependently of the adiabatic component, and it sources the curvature perturbation. In the next 
section, we will derive the corresponding third-order equations, which are needed for the study of 
the primordial trispectrum of curvature perturbation. 
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,3 



(5p^(64) 



III. DERIVING THE THIRD ORDER EQUATIONS OF MOTION 

Let us start by identifying the third order curvature perturbation on a hypersurface of uniform 
energy density. By expanding Eq. (19), one finds, after some manipulations, that one can write 

with 

p p V p J \ p ^p^ J p^ \ p J 

Moreover, by applying the gauge transformation at third order [38] 

<5T(3) ^ <5T(3) + C^^^^ST^^) + (i^l^^ + /:e(.))^T(i) + {C^^^^ + C^^^^C^^^^ + ^4,^)T(o) (65) 

to the covector C,a (for which the formula simplifies since = 0) , one verifies that C^^) indeed 
gauge invariant at third order on large scales. From its definition Eq. (64), one sees that it has the 
desired interpretation as the third-order curvature perturbation (or perturbation in the integrated 
expansion a) on uniform energy density hyper surf aces. It also matches the definition of the same 
quantity on large scales given in [39]. Its evolution equation can be found by expanding Eq. (20): 

C(3)' = _^r(3) - 3^r(2)cW' - 3^r(i)(cw? - ^r^k^^ (66) 

where F^^) is defined in the same way as C^^) gqg (g3) ^nd (64) with the replacements a —>■ 

As the non-adiabatic pressure perturbation F^ is ultimately related to the presence of the 
entropic perturbation (see Eq. (38)), one needs to identify the latter at third order. This is similar 
to the analysis carried out for the curvature perturbation, except that has to bear in mind that 
6a^^^ does not transform as the second-order perturbation of a scalar quantity. Instead 

<5a(^) ^ + e^,f' + 1^0,) (e^.^a')' + ^^^Sa' - f^./Ss (67) 

where the last term is different from the transformation properties of say 6p^'^\ Then, using 
Eqs. (117) and (119), one finds 

a' a' 2a'^ 2a'^ 

+ (^^±^Ss9M, (68) 



2a 
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where Ss^^^ is a gauge invariant quantity at third order on large scales, with the explicit definition 

a' ^ 'a' 2(7 ^ (T Dcr' \ cr / 



-;^5s(5s'(5c7. (69) 



2a'- 

Of course, one could have added additional gauge-invariant terms to the definition of 8s^^ ^ e.g. 
terms involving cubic combinations of the entropy perturbation at first order. Our definition is the 
most natural as it represents the purely third order perturbations of the fields perpendicular to the 
trajectory in the approximate comoving gauge where Sa^^^ = Sa^'^^ = 0. One can also build up Ss^^^ 
by identifying the gauge-invariant quantity that precisely reduces to esj6(f>^^^^ in the approximate 
comoving gauge. Using the transformation rules for the field perturbations, this leads to the same 
definition Eq. (69). 

There is more freedom for the adiabatic perturbation as one is not guided by gauge-invariance 
requirements. Expanding Eq. (35) using Eqs. (116) and (118), one obtains 

<5af ) = 9^5(7(3) + ^ ( SadiSs^^'^ + da^'^'^dids] + ^(6s' + 9'5a)di8s' + -^{5s' + O'SafdiSa 
a' \ J a ^ 2(7 ^ 



with 



5(7(3) ^ e,, V(3) + ^{5s6s(^y + Ss^^^Ss') + J^6s^6s' (71) 

2(7 6(7 ^ 



and 



F.(3) ^ 1 (^5s^^)diSs' + SsdM'^y - Ss^^yOiSs - Ss'dM''^) , (72) 

which is the natural generalization to third order of the non-local term Vi. The global consistency 
of the formalism leads us to the particular choice Eq. (71) for Sa^^^ as will be clear from the 
subsequent discussion of the momentum. 

Expanding Eq. (50), one finds, after a long calculation, 

\ a a 2a 

1 fa"9' e"\ 2 jVss + e'') , 2. ^ , X /X 



2a' 2a ^ \a' J 2a ^ 2a' I 3 I 

5e^'^^di5a SedM^^^ SaSe'djSa Se^diSa a" , , ^ , 

H ^TTT^ h rrj^oeoaOiOa 



a 



2 ' 2a'3 ' a'^ 
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a' a ^ a' 2a' 



^vP + '^6sV, + ^,V, (73) 

OCT (T 

The third-order momentum is similar in form to the second-order one, Eq. (52). In particular, it 
also cannot directly be written as a total gradient. Of particular relevance are the last three terms. 
While the last two involve the second-order quantity V^, negligible in the approximate comoving 

(3) (3) 

gauge at second order, there also appears its third-order generalization . again vanishes 
when the total entropy perturbation factorizes in terms of its time and spatial dependence, and we 
neglect it for the same reason that we neglect Vi at second order, namely that such differences of 
spatial gradients are heavily suppressed on large scales in both inflationary and ekpyrotic models. 
When both Vi and V> are negligible, one can define the comoving gauge at third order on large 
scales as Sa = Jcr^^) = 6a^^^ = 0. 

Using the definitions above, we can define the (approximate) comoving energy density at third 
order. Requiring that this should be a gauge invariant quantity at third order on large scales, 
reducing to the third-order perturbation of the energy density in the approximate comoving gauge, 
we are led to define 

,2^ 



a' a' \ a' J \ a' 2g ^ ) a' \ a' J 



- ^ ( V' - 'i^A - 1 { I)' ^-'^'^^ - ^ ( I)' - h 



1 IP 



ipd' ^ ,^ p'd'^ ^ 3 p' 



5a' 



_ ^(<5s(5a(2) + 5s(2)<5a) - ^5p'5s5a + tt^5a5a'5s - ^Sa^ -^SaSs'^ 
a ^_ _ <Z _ ^1 2(7 3 _ 

+ ^373 5s 5a + + ^ - J '^^^^ - -^^^^^ ^'^ - -^^^^^^s. (74) 

Expanding Eq. (37), with the experience of the calculation of one obtains, using Eqs. (120), 
(121), (122) and (123) 

,(3) = a,6(3) + ^5,6(2)' + ^5,6(^)' - -^5a'die('y + ^9,6(1)" 



5a 
+ — 



2a'3 2a'2 

3V, ( la" + 3hA + - ^ i5ed,5s - 5sd,5e) 

I a a 



3HV^''-'J$SsV. + y^ (75) 



'9 

a ^ 



where the two non-local terms Vi and V^*-^-* naturally appear. 



We now have all the necessary ingredients to derive the third-order equations of motion on large 
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scales. Using the first and second-order results (which can be read off from Eqs. (59) and (61)) 

r(^) « -2Vs5s, (76) 
r(2) « -2y,,5s(2) - VssSs'' + ^^SsSs' , (77) 

we deduce from the equivalent of Eq. (63) that 

r(3) ^ -2V;,5s(3) _ 2V,,5s5s^^^ + Y^{5s^^)5s)' - U^^sss^s^ + ^SL^^^^Ss' 5s' . (78) 

Hence, from Eq. (79), we obtain 

+^ (^,^5.5.(2) - ^,VsV,.5sHs' + \VsVss5s^) + ^^Vps' . (79) 

This equation describes the evolution of the curvature perturbation on large scales at third order, 
and how it is sourced by the entropy perturbation up to third order. It is the third-order counterpart 
of Eqs. (59) and (61). Note that, once the potential V becomes irrelevant, C is conserved on large 
scales. This is important, as the potential indeed becomes irrelevant in the approach to the big 
crunch in ekpyrotic models. 

To complete the system of equations, one needs to derive the second order (in time) equation 
of motion for Ss^^\ This is done by expanding Eq. (39) up to third order, with the result that 

+ (^2^- + 2^ - 3H^-^ {6s('Usy + (v,,ss - 10^ - 18^ j Ss^'^Ss 
+^Ss'' + f % + 3^ + 3H^ - 2% - 6^] 5s'Hs 

/ e'e" 3 - v„Vss o"'v^ Vss o"'\ , , 

+ -10— - —Vssa - 5^^ - l^jf - + UH^ ds'ds^ 

\ a ^ 2(7 a ^ a a ^ ^ } 

+ ^Vssss - ^-Vsss + 2^ + 21^^ + 27^72 M^' = ■ 80 
\6 6 a' a a a ^ j 

Similarly to the second-order calculation, one can derive the same equation by neglecting the 
covector Ca in Eq. (39) from the beginning, if at the same time one assumes that both = and 
V^'^ = 0. 

Eqs. (79) and (80) are the key results of this section. They describe the coupled third order 
perturbations of a system of two scalar fields on large scales without any approximations. In the 
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following sections, we will use them to evaluate the primordial trispectrum in various ekpyrotic 
models of the universe. However, we emphasize that these equations are valid on superhorizon 
scales whatever the dynamics of the scale factor of the universe, so that they are equally well 
suited to the study of cosmological perturbations generated during a period of inflation. 

IV. NON-LINEARITY PARAMETERS AT THIRD ORDER 

The observable that is relevant for comparison with observations is the curvature perturbation 
^, and hence this is the variable that we will focus on. In particular, we wish to evaluate the 
trispectrum, i.e. the Fourier transform of the 4-point function of the curvature perturbation. We 
are interested in calculating the classical non-linearity generated on super-Hubble scales. This was 
shown to be a good approximation for the bispectrum in [19], and consequently we expect this to 
be a good approximation for the trispectrum also. In other words, we expect the non-linearity on 
sub-Hubble scales to be subdominant compared to the non-linear classical super-Hubble evolution. 

We only need to consider the connected part of the 4-point function, i.e. the part that is not 
captured by products of 2-point functions. If a distribution is exactly gaussian, all information is 
contained in the 2-point function. In particular this means that when n is odd, all n-point functions 
are zero, and when n is even they can be decomposed into products of 2-point functions. Here we 
therefore consider the part of the 4-point function which is not already captured by the product of 
two 2-point functions. Thus we define 

(^1^2^3^4)0 = i^TrfS^ki + k2 + k3 + k4)T(fei, ^2, h, h), (81) 

where T describes two different shape functions parameterized by the A;-independent parameters 
tnl and gNL, see [8] for more details. These are defined via 

54 

T = TArL[-P(fci3)-P(fc3)-P(fc4) + llpermutations]-F— 5iVL[-P(/c2)^'(A;3)P(A;4) + 3permutations], (82) 
where P{k) is the power spectrum defined by the 2-point function 

(CfciCfe) = {27rfS\ki + k2)P{h), (83) 

and kjj = kj + kj . 

The parameters rjvi and (/jv^ can be related to an expansion of the curvature perturbation in 
terms of its gaussian component, as follows. The Bardeen space-space metric perturbation [40] 
is expanded as 
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where denotes the Hnear, gaussian part (up to an unimportant zero mode such that = 0). 
During the era of matter domination, this variable is related to the curvature perturbation via 
C = f'J'-ff! so that we obtain the expansion 

C = CL + lfNLCl + -^gNLCl (85) 

with (^L being the linear, gaussian part of i^. tjvl is directly related to the square of /atl, explicitly 

36 

TNL = ^fNL- (86) 

The strategy for calculating the non-linearity parameters is straightforward: first we solve the 
equations of motion for the entropy perturbation up to third order in perturbation theory. This 

allows us to integrate the equation of motion for Q, also at the first three orders in perturbation 
theory, and then we obtain the non-linearity parameters by evaluating 

5 /t*'C^'^' 

25 f/'C^'^'' 

^^■^ = "9"(J7^(Iyji' ^^^^ 

where the integrals are evaluated from the time tj that the ekpyrotie phase begins until the eon- 
version phase has ended at and C, has evolved to a constant value. The next sections will deal 
with the details of these calculations. 



V. EKPYROTIC EXAMPLES 



In ekpyrotie and cyclic models of the universe, the standard cosmological puzzles, such as the 
horizon and flatness problems, are resolved by a phase of slow contraction before the big bang [41]. 
During this ekpyrotie phase, the equation of state Wek (the ratio of pressure to energy density) is 
very large, with the consequence that the ekpyrotie matter (usually modeled as a scalar field with a 
steep and negative potential) quickly comes to dominate the energy budget in the universe. In other 
words, the contributions to the total energy density due to spatial curvature and anisotropics, for 
example, become suppressed, and in this way the flatness problem is solved. The horizon problem 
does not even arise, as there was plenty of time before the big bang for different regions of the 
universe to be in causal contact with each other. 

Interestingly, and just as in inflation, the ekpyrotie phase can (quantum-mechanically) produce a 
nearly scale-invariant spectrum of density perturbations at the same time as (classically) flattening 
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the universe. However, early studies [42-44] of single-field ekpyrotic models uncovered a subtlety 
related to the fact that the perturbations are generated during a contracting phase: namely, the 
mode that picks up a scale-invariant spectrum corresponds to a small time-delay perturbation to 
the big crunch, and this mode turns out to be the decaying mode in a contracting universe. The 
adiabatic mode, on the other hand, picks up a blue spectrum. Hence, unless these two modes mix 
during the big crunch - big bang transition (which they might do due to higher-dimensional effects 
for example [45, 46]), the spectrum of the growing mode during the subsequent expanding phase 
will be incompatible with observations. 

However, it is unnatural to restrict the analysis to a single scalar field. In models inspired by 
heterotic M-theory for example [47-50], there are two universal scalars, namely the radion mode 
and the volume modulus of the internal Calabi-Yau manifold. The radion determines the distance 
between two orbifold planes, and in this braneworld scenario, our currently observed universe is 
identified with one of the orbifold planes. In this picture, the ekpyrotic potential corresponds to 
an attractive force between the orbifold planes, which eventually causes the branes to collide. This 
collision then looks like a big bang to an observer on the brane. Quickly afterwards, the distance 
between the branes becomes fixed again, but the branes themselves expand during the phases of 
radiation, matter and dark energy domination. Actually, the dark energy can be modeled by the 
same inter-brane force that is at play during the ekpyrotic phase, it just corresponds to a rather 
flat plateau at large field values. Inexorably though, the branes are attracted to each other again, 
a new ekpyrotic phase takes place, which is then swiftly followed by a new brane collision, etc. In 
this way a cyclic picture naturally emerges [51]. 

But back to the two scalar fields: the theory we are considering can be described by the effective 
action [52] 

s = J V^[R-l{d<Pif - ^{dh? - v{cp,,M, (89) 

where (f)i and ^2 are related by a field redefinition to the radion and the Calabi-Yau volume 
modulus. We are assuming that during the ekpyrotic phase, both fields feel an ekpyrotic-type 
potential, e.g. 

V{<f>i,<f>2) = -Vie-^^^' - Fse-"^'^^ (90) 

In fact, as discussed in detail in [53], it is much more natural to discuss the evolution in terms 
of the variables a and s pointing transverse and perpendicular to the field velocity respectively 
(these are the variables we have used in sections II and III). In terms of those variables, the above 
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potential can be re-expressed as an ekpyrotic potential for a combined with a tachyonic potential 
in the transverse direction. The evolution proceeds along this unstable ridge in the potential. This 
instabihty is necessary for the amplification of scale- invariant entropy perturbations [54] , and leads 
to an interesting global structure [5]. 

Recasting the potential in terms of a and s in fact suggests that, in the absence of a microphysical 
derivation, we should use a more general, parameterized, form for the potential instead of (90). 
The feature that we want to keep is that we should still obtain a nearly scale-invariant spectrum 
of entropy perturbations. Hence, during the ekpyrotic phase, we require the potential to have the 
following features: 1. Along the a direction, it must be of ekpyrotic form, i.e. negative and steep; 
more precisely, we want V^a ~ — 2/f^. 2. Along the s direction, it must be tachyonic, i.e. we need 

= and Vss < 0. As we just said, this instability is necessary in order for entropy perturbations 
to grow. 3. In order to obtain a nearly scale-invariant spectrum of entropy fluctuations, we must 
have Vaa ^ ^ss- Small deviations from these relations lead to small differences for the tilt of the 
spectrum of fluctuations [53]. Thus we write the potential as 

Ve, = -Voe^-[l + es' + ^e'/'s' + ^eh^ + •••], (91) 

where we expect ks,K4 ~ 0(1) and where e ~ O(IO^) is related to the equation of state Wek 
via e = 3(1 -|- W(,k)/2. For exact exponentials of the form (90), one has K3 = 2\/2(cf — c^)/\c\C2\ 
and = 4:{c\ + C2)/(cfc2(cf + c^)). Before continuing, we should comment on a subtlety: strictly 
speaking, s is not a field, since by definition s = along the background trajectory. Hence, when 
we expand the potential in terms of s, we simply use s to denote the direction transverse to the 
ridge of the potential. We can only use the more compact notation here because the background 
trajectory is a straight line during the ekpyrotic phase, and hence the meaning of s is unambiguous: 
in terms of the original fields, one could use the replacement s = (c202 — ci4>i)/ \/ cf + c|, and this 
would give identical results. 

Along the ridge, the background scaling solution is given by 

a{t) = {-tf' a = -^lii [-^/eVot^ s = 0, (92) 

where time runs from large negative values at the beginning of the ekpyrotic phase to i = at 
the big crunch. This scaling solution corresponds to a straight line trajectory in scalar field space. 
At the same time, nearly scale-invariant entropy perturbations are produced and amplified, their 
final amplitude being determined by the depth of the potential. Subsequently, any bending of 
the trajectory in field space will result in the conversion of entropy perturbations into curvature 
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(adiabatic) perturbations, which inherit the same spectral index [13]; accordingly, this mechanism 
is termed the entropic mechanism. 

A bending of the trajectory is expected in many models on generic grounds - we will discuss two 
cases in detail in sections VC and VD. The first is where the bending occurs during the kinetic 
energy dominated phase that leads up to the big crunch ("kinetic conversion"), and the second is 
where the trajectory starts out slightly off the ridge, so that after some time, but still during the 
ekpyrotic phase, it naturally rolls off the ridge to one side ("ekpyrotic conversion"). As we will 
see, both the sign and the amplitude of the non-linearity parameters depend rather crucially on 
the details of this conversion process. 

Before continuing, we should state an important assumption that we are making in this work: 
namely, we assume that the curvature perturbation passes through the big crunch - big bang tran- 
sition essentially unchanged. The reason for doing so is that the perturbations we are considering 
are vastly larger than the horizon size around the time of the crunch, and hence, due to causality, 
it seems reasonable to assume that long-wavelength modes suffer no change. Nevertheless, in the 
absence of a complete quantum-gravitational treatment of the big crunch - big bang transition, 
this statement remains an assumption subject to possible revision in the future. 



A. Initial conditions for the entropy mode 

During the ekpyrotic phase, the trajectory is a straight line with 9' = 0, (see Eqs. (30) and (32)) 
and so the equation of motion for ds^^^ simplifies to 

+^Ss'' + f % + 3^ + - 2%^ Ss'^Ss 



+ ( -^,Vssa - - 3i/|f ) 5s'5s' + (^-V,ssss + 2^ J 5s' = Q. (93) 

Using the following useful expressions, valid during the ekpyrotic phase, 

- _ V2 - _2 - _ - _ 2V2€ 

y - V --- V --^ V --^ 

v,s — o v^ss — ^2 ~ -f-2 ^,ssss — ^2 ' 
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together with the first- and second-order entropic equations (60), (62) in the hmit when 0' = 0, we 
iteratively find that the entropy perturbation (to leading order in 1/e) can be expanded as 

Ss = 6sL + '^6sl + e(S + S - hssi, (94) 
8 ^60 80 S'^ V y 

with 5s L oc 1/t. The above equation specifies the initial conditions for the start of the conversion 
phase. 



B. Integrated contributions during the ekpyrotic phase 

There are contributions to both /jvl and g^L even during the ekpyrotic phase, since some of 
the terms in the equation for ^' do not vanish when 9' = 0. Specifically, during the ekpyrotic phase, 
we have (C here is the full C^^^ + C^^^ + C^^^) 

a ^ a' 3 



— 5s 5s5s' — 5s 

t 2 6t 

5s1 llKs^/e5sj^ 

2i m 



(95) 



This leads to 



_ 5 [5sL{tek-end)? .q^s 

_ 275 A^SL{tek-end)f /^'7^ 
9N Lintegrated — j^gg'^^Ve [(^{tf)]^ ' 

where tf denotes the time at the end of the conversion stage, when ( becomes constant on large 
scales, and tek-end denotes the time at the end of the ekpyrotic phase. 

Numerically speaking, both J n Lintegrated and QMLintegrated are typically subdominant contribu- 
tions to the total non-linearity parameters. 



C. Kinetic conversion 



In the original ekpyrotic and cyclic models, the phase dominated by the steep, ekpyrotic po- 
tential Vgfc comes to an end (at t = tefc-end < 0) before the big crunch/big bang transition (at 
i = 0), and as the ekpyrotic potential becomes negligible the universe becomes dominated by the 
kinetic energy of the scalar fields. In this subsection, we consider the case where the conversion 
from entropic to curvature perturbations occurs during this kinetic energy dominated phase. This 
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FIG. 1: /nl as a function of the duration of conversion and for the values K3 = —5,0,5 and e = 36. In 
each case, we have plotted the results for four different reflection potentials, with the simplest potentials 
((/)^^, (sinh (/)2)^^) indicated by solid lines, while the dashed ((sinh(/)2)^^+(sinh(/)2)^^) and dotted (02^^+(/'2^^) 
lines give an indication of the range of values that can be expected. As the conversions become smoother, the 
predicted range of values narrows, and smooth conversions lead to a natural range of about —50 < /nl ^ +60 
or so. 



occurs naturally in the heterotic M-theory embedding of the cyclic model [55] (briefly discussed 
towards the beginning of this section) because the negative-tension brane bounces off a spacetime 
singularity (see [14] for a detailed discussion) - thus creating a bend in the trajectory in field space 
in the 4-dimensional effective theory [52]- before it collides with the positive-tension brane (the 
latter collision corresponds to the big crunch-big bang transition). This bending of the trajectory 
automatically induces the conversion of entropy to curvature perturbations [13]. From the view- 
point of the 4-dimensional effective theory, the bounce of the negative-tension brane corresponds 
to the trajectory in scalar field space being reflected off a boundary of scalar field space, and it can 
be modeled by having a repulsive potential Vrep in the vicinity of the boundary. Here we take the 
boundary to be located at (f)2 = 0, and we only consider conversions during which 6' > ~ other 
cases can be related to these by an appropriate change of coordinates. The repulsive potential can 
in principle be calculated, given a specific matter configuration on the negative-tension brane [14]. 
Here, in order to be general, we consider four different forms for the repulsive potential that we 
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-25000 ' ' 

FIG. 2: gjvL as a function of the duration of conversion, with K3 = K4 = and for four different reflection 
potentials, with the same line style assignments as in Fig. 1. As the conversions become smoother, the 
predicted range of values narrows considerably, allowing us to make rather definite predictions. 

consider to be representative, namely 

Vrep oc 02 ^ 02 ^ + </'2 ^ (smh02)~^ (slnh 02)"^ + (sinh02)~^ (98) 

with the overall magnitude adjusted in order to obtain various values for the duration of the 
conversion (see below) . These potential forms should give an indication of the range of values that 
one can expect the non-linearity parameters to take. 

An important parameter turns out to be the duration of the conversion, measured by the 
number of Hubble times during which most, say 90 percent, of the conversion takes place, i.e. the 
duration is measured by the number of e-folds by which the scale factor shrinks during conversion. 
Conversions lasting less than about 0.2 Hubble times correspond to what we call sharp conversions, 
while those that last on the order of 1 Hubble time correspond to smooth conversions. For /tvl, 
the range of predicted values is considerably narrower as the conversion becomes smoother [17, 18], 
as illustrated in Fig. 1. It is possible to estimate Jnl semi-analytically, with the result that for 
smooth conversions [18], 

fNL ~ ^ + 5. (99) 
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FIG. 3: This figure shows g^vi to be proportional to e, i.e. proportional to the equation of state Wek- 



The values for ttvl = 36/25/^^ follow directly from these calculations, and hence they do not 
require an extensive discussion. 

In order to calculate qnl, we have solved and integrated numerically the equations of motion 
(80) and (79), using the expression (94) as the initial condition for the entropy perturbation. C, is 
initially zero. The initial conditions for the background trajectory are taken to be those suggested 
by the heterotic M-theory embedding [52], i.e. we take the incident angle of the trajectory with 
respect to the boundary surface </)2 = (and thus the reflection potentials) to be 7r/6. On changing 
this angle, we expect our results to change only modestly [18]. 

The results of our numerical calculations are shown in Figs. 2 to 5. In each case, we have 
plotted the results obtained for the four repulsive potentials (98). Fig. 2 shows that, even more 
so than for the calculation of Jnl the range of predicted values for gjvL narrows drastically as the 
conversion process becomes smoother. In fact, for sharp conversions, typical values are very large 
in magnitude, and we expect these to be observationally disfavored shortly, if they aren't ruled 
out already. Thus, phenomenologically speaking, it is much more interesting to focus on smooth 
conversions. In this respect, it is important to note that for smooth conversions, the predicted 
range of values for ff^i lies within the current observational bounds [56, 57]. Hence, we will now 
focus on conversions that occur while the scale factor evolves by more than about half an e-fold. 
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FIG. 4: This figure indicates that gjqi, depends hnearly in K4, the parameter we are using to specify the 
fourth derivative of the ekpyrotic potential with respect to the entropic direction. 

The next three figures show the dependence g^L on the potential parameters ^3,^4 and 
e = 3(1 + ti;efc)/2, as defined in (91). Fig. 3 shows that gj^n is proportional to e, and we have 
verified that this result holds independently of the values of K3 and K4. Similarly, Fig. 4 indicates 
that gNL depends linearly on K4 (independently of e and K3), while Fig. 5 indicates that g^L 
depends approximately quadratically on K3 (independently of e and K4). In fact, all our numerical 
results indicate that g]\ji scales with e, K3, K4 exactly as the third-order coefficient in the expression 
(94) for the entropy perturbation during the ekpyrotic phase, and we have found that we can fit 
all our data by the approximate formula 



i.e. gNL is roughly speaking directly proportional to the third-order coefficient in the expansion 
(94) of the entropy perturbation during the ekpyrotic phase. This formula is one of the key results 
of this paper. It shows that gNL is proportional to the equation of state of ekpyrosis w^k ~ £• 



(99), we can deduce that a large I/atlI tends to make gNL more positive due to its dependence on 
K3. Most significant is perhaps the term independent of K3, K4. This term is negative and typically 
of order 1000, and it implies that gNL can be expected to be negative unless \nz\ is large and/or 




(100) 



Hence it scales in the same way as ttvl oc oc (-v/e)^. Also, from the fitting formula for Jnl 
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FIG. 5: QNL can be seen to depend approximately quadratically on K3, the parameter determining the third 
derivative of the ekpyrotic potential with respect to the entropic direction. 



K4 is large and positive. In fact, the smaller \fNL\ is, the more 'tuned' the value of K4 has to be 
in order to make qnl positive - for K3 ~ 0, we would require K4 ~ 25 just in order to make g^L 
zero. Hence, any accidental degeneracy at the level of fNL between simple inflationary models (in 
which both f^i and qnl are expected to be very close to zero) and cyclic models is extremely 
likely to be broken at the level of the trispectrum. Our result indicates that it is far from true 
that "anything goes", e.g. it would be natural to find Jni = 20 and g^L = —1500, while it 
would be very unnatural to get Jnl = 20 together with gi\ji = +3000. Thus, for the currently 
best-motivated cyclic models in which the perturbations are converted during the kinetic phase, 
Eq. (100) provides a distinct observational imprint. 

For completeness, we should note that there is the possibility that the dominant amount of 
conversion could occur after the big bang, via the mechanism of modulated preheating suggested 
by Battefeld [58] . The idea here is that massive matter fields could be produced in large quantities at 
the brane collision and could dominate the energy density right after the big bang. If these massive 
fields were to couple to ordinary matter via a function h{5s), then their decay into ordinary matter 
would occur at slightly different times for different values of 6s. In this way, the ordinary matter 
perturbations would inherit the entropic perturbation spectrum. This conversion would typically 
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also happen during a phase of kinetic energy domination, so that we would expect the orders of 
magnitude of /jvl and g^L to be similar to the examples above. However, in the absence of a 
detailed model, which would determine h{Ss), we cannot do better than this order-of-magnitude 
estimate. 

D. Ekpyrotic conversion 

In this section, we analyze the case where the conversion of entropy into curvature modes 
occurs during the ekpyrotic phase. Indeed, if the background trajectory starts out slightly off the 
ridge in the ekpyrotic potential, then after some time, the field will naturally roll off down one 
side [15, 59, 60]. For this scenario, the non-gaussianity up to second order has been analyzed by 
Koyama et al. [19] using the SN formalism, with the result that /jvl = — 5/12cf, for the case 
where 01 is the field that is frozen in at late times, and where ci is the constant appearing in 
the ekpyrotic potential (90). For ekpyrotic conversion, the calculation is most easily performed, 
and the result most easily expressed, in terms of the potential (90), which is why we are adopting 
this restricted form here. In working with a parameterized potential like (91), the bending of the 
trajectory can be more complicated, in the sense that there can be multiple turns, and one has 
to decide when to stop the evolution. In this case, the results are strongly cut-off dependent, and 
without a precisely defined model specifying the subsequent evolution, it is impossible to make 
any generic predictions. For kinetic conversion, this problem does not arise, since the refiection 
potential is entirely unrelated to the potential during the ekpyrotic phase. The bottom line is that 
we will stick with the form (90) in this section. 

The 5N formalism is particularly well suited to the case of ekpyrotic conversion, as the back- 
ground evolution is simple. In fact, it turns out that by making the approximation that the bending 
is instantaneous, it is very easy to find an approximate formula for the non-linearity parameters 
at any chosen order in perturbation theory. Below, we will summarize the calculation of [19] and 
extend it to third order. Subsequently, we will compare the approximate formula thus obtained 
with the result from solving and integrating the equations of motion numerically. In a sense, this 
provides a check of the instantaneous bending approximation. At second order in perturbation 
theory, it has been checked explicitly that if the instantaneous bending approximation is relaxed, 
then both methods (namely a numerical calculation using the 5N formalism, and a numerical cal- 
culation using the equations of motion directly) agree to high precision in the value of /nl that 
they predict, and moreover that the results are in good agreement with the formula quoted above 
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[18]. 

In order to implement the SN formalism, we have to calculate the integrated expansion N = 
J Hdt along the background trajectory. Initially, the trajectory is close to the scaling solution 
(92). Then, we assume that at a fixed field value Ssb away from the ridge, the trajectory instantly 
changes course and rolls off along the ^2 direction. At this point, the trajectory follows the single- 
field evolution 

a{t) = (-t)"^^"^ (j)2 = — In(-t) + constant 0i = constant. (101) 

C2 

It is now easy to evaluate the integrated expansion, with the result that 

2 

N = — ^ In + constant, (102) 

where Hb denotes the Hubble parameter at the instant that the bending occurs. Note that all 
C2-dependence has canceled out of the formula above. At the end of the conversion process, we are 
interested in evaluating the curvature perturbation on a surface of constant energy density. But, in 
comoving gauge, the curvature perturbation is equal to a perturbation in the integrated expansion. 
Then, if we assume that the integrated expansion depends on a single variable a, we can write 

( = 6N = Nc^Sa + ^Naa{5af + ^N^aU^af- (103) 

In our example, we indeed expect a change in N to depend solely on a change in the initial value 
of the entropy perturbation 5s. Now, from Eq. (94), we know that 6s l oc 1/t oc H, and hence we 
can parameterize diflFerent initial values of the entropy perturbation by writing 

5sL = aH. (104) 

Note that since Ssl is gaussian, so is a. With this identification, we have 

5s = aH+ '^{aHf + e(^ + M _ ^)(ai^)^ (105) 
8 ^ ^ ^60 80 5^^ ^ ^ 



so that at the fixed value (Js = we have 



Now we can immediately evaluate 



1 . X 

a oc — . (106) 
Hb 
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Cl 


C2 


fNL,SN 


tnl.sn 


gWL.SN 


Inl 


tnl 


9nl 


10 


10 


-41.67 


2500 


2315 


-39.95 


2298 


2591 


10 


15 


-41.67 


2500 


2315 


-40.45 


2356 


2813 


10 


20 


-41.67 


2500 


2315 


-40.62 


2377 


3030 


15 


10 


-93.75 


12660 


11720 


-91.01 


11930 


13100 


15 


15 


-93.75 


12660 


11720 


-92.11 


12220 


13830 


15 


20 


-93.75 


12660 


11720 


-92.49 


12320 


14440 


20 


10 


-166.7 


40000 


37040 


-162.5 


38020 


41320 


20 


15 


-166.7 


40000 


37040 


-164.4 


38930 


43170 


20 


20 


-166.7 


40000 


37040 


-165.1 


39240 


44490 



TABLE I: Ekpyrotic conversion: the values of the non-linearity parameters estimated by the 5N formalism 
compared to the numerical results obtained by directly integrating the equations of motion. 



and, similarly 



.2^2 



N 



.2^3 



In this way, with very little work, we can estimate the non-Hnearity parameters 



Tnl 

gNL 



5N, 
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12 



25N 



,cxcxcx 



5m% 



25 

los' 



(108) 

(109) 

(110) 
(111) 



We are now in a position to compare these estimates to the numerical results obtained by solving 
and integrating the equations of motion (80) and (79). In order to do this, we choose initial 
conditions that are given by the scaling solution (92), except that we increase the initial field 
velocity \(f)2\ by 0.1 percent. This causes the trajectory to eventually roll off in the ^2 direction, 
and to quickly approach the single-field solution (101). The results for several values of ci and C2 
are shown table I, alongside the values estimated by the SN formulae. 

It is immediately apparent that the general trend is accurately captured by the 6N formulae. 
However, one may notice that the agreement is slightly less good at third order than at second, 
and also, that the 6N formulae tend to slightly over-estimate ttvl and slightly under-estimate gNL- 
But given the quickness of the SN calculation and the complexity of the third order equations of 
motion, the agreement is pretty impressive. Of course, the SN formula was derived subject to the 
instantaneous bending approximation. Without this approximation, we would expect a numerical 
scheme that uses the SN formalism to yield results in close agreement with our numerical results. 
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The general conclusion is that, contrary to /jvl, the sign oi qnl turns out to be always positive. 
Moreover, both tjvl and g^L scale very fast with increasing equation of state Wek ~ cf , and hence 
we can expect future observations to be highly constraining for this type of conversion (current 
observations are in fact already rather constraining regarding /nl [56, 57]). 

Finally, we should note that a conversion mechanism similar to ekpyrotic conversion (in the 
sense that the conversion is assumed to happen during the ekpyrotic phase) has been considered 
in the new ekpyrotic models of Buchbinder et al. [16, 53, 61], except that the roll-off from the 
ridge is expected to be caused by a feature in the potential rather than by initial conditions that 
are slightly off-centered. We would simply like to note that for these models, there is in fact 
considerably more flexibility, as the ekpyrotic phase is followed by a ghost condensate phase. The 
details of how the transition to the ghost condensate phase occurs will determine what one expects 
for the non-linearity parameters fNL,TNL,gNL- For example, it is conceivable that the ekpyrotic 
phase might be followed by a kinetic phase during which the ghost condensate starts dominating. 
In that case, one would expect the results to be more closely aligned with those presented in the 
previous section. In the absence of a concrete model, we must postpone making definite predictions 
for these models. 

VI. DISCUSSION 

The analysis of the distribution of density fluctuations, either via the CMB or via large-scale 
structure surveys, currently offers the best prospects for gaining information about the time around 
the big bang. In ekpyrotic and cyclic models, the pattern of these fluctuations is imprinted during 
a slowly contracting ekpyrotic phase preceding the current expanding phase. We have restricted 
our analysis to the most robust mechanism to date by which the density fluctuations can arise in 
these models, namely the entropic mechanism in which entropy perturbations are generated first, 
and are then converted into curvature perturbations just before the big bang. The conversion can 
happen in at least two distinct ways, either directly during the ekpyrotic phase or during the kinetic 
energy dominated phase that leads up to the big crunch - big bang transition. In this paper, we 
have assumed that the dominant amount of conversion occurs before the big crunch, and that the 
resulting density perturbations evolve essentially unchanged up to the time of nucleosynthesis, so 
as to become the 'primordial' density perturbations. 

For both conversion modes, we have found that the combined consideration of the non-linearity 
parameters /jvl (and by extension tnl) and g^L results in a distinct observational imprint, that 
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should enable one to select or rule out these models on observational grounds in the foreseeable 
future. In short, for ekpyrotic conversion, /^vl is always negative and typically of 0(10 — 100), 
while gNL is always positive and typically of 0(10^ — 10^). For kinetic conversion, on the other 
hand, /tvl can have either sign and is typically of 0(10), while qnl is typically negative and 
of 0(1000). Thus, typical values arc substantially different from the predictions of simple single- 
field inflationary models, and the local form of the produced non-gaussianity should also easily 
distinguish ekpyrotic models from single-field inflationary models with non-canonical kinetic terms. 
The comparison with multi-field inflationary models is more subtle, as some of those models can 
allow for virtually any values of /atl and gNL- However, should the observed values happen to 
lie in the ranges predicted in this paper, then ekpyrotic models will provide strong candidates for 
a model of the early universe due to the correlation between the values of /jvl and gNL- It will 
certainly be exciting to compare the predictions derived here to the results of future observations. 
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Appendix - Useful formulae 



yP = ^=1- - + 1^(1)2 - ^(3) + 3^(2)^(1) - (112) 

OT 2 2 

= 0, (113) 

where A is the lapse function defined by ^foo = — (1 + 2A). 



Cf^ = C!f^' -Adf^' -^C!P' + \A^C:f^' (114) 



5ei = j^{5s' + e'5a)el (116) 
5el=~{5s' + 9'5a)ei (117) 
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2a 



5a 



+ ^,(2)' + ^ 5,' + Ua 



a' 



a' 



- {pa' - 9' 6s) {Ss' + e'Sa) + 6' (^Sa^'^'^ - ^Ss6s'^ 



(118) 



= {6s' + e'Saf -^f ~ {^c^' - 0'6s) {5s' + 9' 5a) + 9' - -^5s5s'^ 



6a 



+ Ss^^y + r4r 6s' + Ua 



(119) 



<5c7(^) = + 5a' - e'5s 



(120) 



um'6s - e'6(^h + {Vss + SP)^^ - ^5s5s' + ^5s5e (121) 

2a 2a 



(1) 



a' a 



a 



(122) 



(2) 



= ^ 

a' 

~ -'2 



o(2)' 



+ 



- A.6a('y 



a' a'2 



_ (v,s + 0'^)— + ^6s6s' - ^6s6e — + -^9' 6s 

2a' 2a ^ a ^ a' p 



(123) 
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In comoving gauge and on large scales 

5Vs « VssSs - ^6s' (124) 

5Vf^ « -YsssSs" + - -^VsaSs5s' - ^Ss'"" - + ^^5s5s') (125) 

,s 2 2a-' ' 2a 2 a' 2a' 

SVss ~ V^,,,5s - 2%<5s' (126) 
a 

_2%(<55(2)' + l-<55(55') (127) 
a 20- 



+^ {-^'Vsss + 3Vsse' + 3^'^) Ss^ (128) 



A fs -2^7(^5 (129) 
« -2^,5s(2) + l-(Vss + Qe'^)Ss^ - ^SsSs' (130) 

Finally, the constraint equations derived in [23] imply, in the comoving gauge and on large scales 

e(^) « 0, (131) 
e(^) ^ 0. (132) 
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